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1 Introduction 

Surely one of the most important and general results of statistieal meehanies 
is the existenee of a relation between the spontaneous fluetuations and the re- 
sponse to external fields of physical observables. This result has applications 
both in equilibrium statistieal mechanics, where it is used to relate the correla- 
tion functions to macroscopically measurable quantities such as specific heats, 
susceptibilities and compressibilities, and in nonequilibrium systems, where it 
offers the possibility of studying the response to time-dependent external fields, 
by analyzing time-dependent correlations. 

The idea of relating the amplitude of the dissipation to that of the fluc- 
tuations dates back to Einstein's work on Brownian motion [l]. Later, On- 
sager [3J [3] put forward his regression hypothesis stating that the relaxation of 
a macroscopic nonequilibrium perturbation follows the same laws which govern 
the dynamics of fluctuations in equilibrium systems. This principle is the ba- 
sis of the fluctuation-dissipation relation (FDR) of Callen and Welton, and of 
Kubo's theory of time-dependent correlation functions [4]. This result repre- 
sents a fundamental tool in nonequilibrium statistical mechanics since it allows 
one to predict the average response to external perturbations, without applying 
any perturbation. 

Although the FDR theory was originally obtained for Hamiltonian systems 
near thermodynamic equilibrium, it has been realized that a generalized FDR 
holds for a vast class of systems. A renewed interest toward the theory of fluc- 
tuations has been motivated by the study of nonequilibrium processes. In 1993 
Evans, Cohen and Morriss [5j considered the fluctuations of the entropy produc- 
tion rate in a shearing fluid, and proposed the so-called Fluctuation Relation 
(FR). Such relation was then developed by Evans and Searles [5] and, in differ- 
ent conditions, by Gallavotti and Cohen [7]. Noteworthy applications of those 
theories are the Crooks and Jarzynski relations [Hi iSj . 

In the recent years there has been a great interest in the study of the relation 
between responses and correlation functions, also with the aim to understand 
the features of aging an glassy systems. While, on one hand, a failure of the 
validity of the equilibrium statistical mechanics (e.g. lack of the ergodicity) 
implies a violation of the fluctuation- dissipation theorem, on the other hand, 
as we will see in Sections |3l [S] [H the opposite is not true: a failure of the 
fluctuation- dissipation theorem (at least in its "naive version" ) does not imply 
the presence of glassy phenomena. It may, for instance, resides in the presence 
of stationary non-equilibrium currents. 

1.1 The relevance of fluctuations: few historical comments 

Let us spend few words on the conceptual relevance of the fluctuations in sta- 
tistical physics. Even Boltzmann and Gibbs, who already knew the expression 
for the mean-square energy fluctuation, 

{{E - {E)f) ^ keT^Cv (1) 



pointed out that fluctuations were too small to be observed in macroscopic sys- 
tems: In the molecular theory we assume that the laws of the phenomena found 
in nature do not essentially deviate from the limits that they would approach in 
the case of an infinite number of infinitely small molecules. . . \10f . . . [the fluc- 
tuations] would be in general vanishing quantities, since such experience would 
not be wide enough to embrace the more considerable divergences from the mean 
values [IT]. 

On the contrary Einstein and Smoluchowski attributed a central role to the 
fluctuations. For instance Einstein understood very well that the well known 
equation ((!]) ... would yield an exact determination of the universal constant 
(i.e. the Avogadro number), if it were possible to determine the average of 
the square of the energy fluctuation of the system; this is however not possible 
according our present knowledge |12j . 

It is well known that the study of the fluctuations in the Brownian motion al- 
lows one for the determination of the Avogadro number (i.e. a quantity at micro- 
scopic level) from experimentally accessible macroscopic quantities. Therefore 
we have a non ambiguous link between the microscopic and macroscopic levels; 
this is well clear from the celebrated Einstein relation which gives the diffusion 
coefficient D in terms of macroscopic variables and the Avogadro number: 

i, . lin, Ml^iiM . ^Z_ , (2) 

t^oo 2t 6NATTr]a ^ ' 

where T and r\ are the temperature and dynamic viscosity of the fluid respec- 
tively, a the radius of the colloidal particle R = NaI^b is the perfect gas 
constant and N a is the Avogadro number. 

The theoretical work by Einstein [1] and the experiments by Perrin [T3] 
gave a clear and conclusive evidence of the relationship between the diffusion 
coefficient (which is measurable at the macroscopic level) and the Avogadro 
number (which is related to the microscopic description). Therefore, already at 
equilibrium, small-scale fluctuations do matter in the statistical description of 
a certain system. 

The main purpose of the present contribution is at showing how such fluctu- 
ations must be accounted for when the system is out of equilibrium. A relevant 
outcome will be the evidence that in order to properly describe small scale out- 
of-equilibrium fluctuations some correlations not present at equilibrium must 
be taken into account. 

The celebrated paper of Einstein, apart from its historical relevance for the 
atomic hypothesis and the development of the modern theory of stochastic pro- 
cesses, shows the first example of FDR. Let us write the Langevin [Tl] equation 
for the colloidal particle of mass m: 



dt V m 

where 7 = QNA'^'na/m, and ^ is a white noise, i.e. a Gaussian stochastic process 
with (C(i)) = and {C{t)C{t')) ~ ^{t ^ t')- An easy computation gives D = 



{V'^)t. Consider now a (small) perturbating force f{t) = FQ{t), where Q{t) is 
the Heaviside step function. It is simple to determine the average response of 
the velocity after a long time (i.e. the drift) to such a perturbation: {6V) = 
F/7. Defining the mobility fi as: {5V) = fiF one easily obtains the celebrated 
Einstein relation (EFDR) 

which gives a link between the diffusion coefficient (a property of the unper- 
turbed system) and the mobility which measures how the system reacts to a 
small perturbation. 

An intcrsting point addressed here is that out-of-equilibrium conditions can 
be hardly obtained from a single Langevin equation when the noise and drag 
terms are local in time. Preserving this locality, and hence the Markovian 
nature of the model, we will see that out of equilibrium a proper description of 
the relation between correlations and responses can be given within an enriched 
space of variables, where also some local field coupled to the variable of interest 
must be considered. This new field is itself a fluctuating object whose dynamics 
is ruled by another stochastic equation. 

2 Generalized Fluctuation-Dissipation relations 

The Fluctuation-Response theory was originally developed in the context of 
equilibrium statistical mechanics of Hamiltonian systems. Sometimes this in- 
duced confusion, e.g. one can find misleading claims on the limited validity of 
the FDR: we can cite an important review on turbulence containg the following 
(wrong) conclusion This absence of correlation between fluctuations and relax- 
ation is reflected in the non existence of a fluctuation- dissipation theorem for 
turbulence |15] . 

On the contrary, as we will discuss in the following, a generalized FDR holds 
under rather general hypotheses [16l [171 HI] ■ 



2.1 Chaos and the FDR: van Kampen's objection 

It has been argued by van Kampen [19] that the usual derivation of the FDR that 
relies on a first-order truncation of the time-dependent perturbation theory, for 
the evolution of probability density, can been severely criticized. In a nutshell, 
using the dynamical systems terminology, the van Kampen's argument is as 
follows. Given a perturbation (5x(0) on the state of the system at time t = 0, 
one can write a Taylor expansion for (5x(<), the difference between the perturbed 
trajectory and the unperturbed one: 

^^^W = E £||'^^^(0) + O(|<5x(0)n . (5) 



Averaging over the initial condition one has the mean response function: 



In the case of the equilibrium statistical mechanics we have p(x) oc exp( -~[3H{x.) 
so that , after an integration by parts, one obtains 






^..(0=/3 -.W^^S^)> (7) 



which is nothing but the differential form of the usual FDR. 

In the presence of chaos the terms dxi{t)/dxj{0) grow exponentially as e'^*, 
where A is the Lyapunov exponent, therefore it is not possible to use the linear 
expansion ([5]) for a time larger than (1/A) ln(i/|(5x(0)|), where L is the typical 
fluctuation of the variable x. On account of that estimate, the linear response 
theory is expected to be valid only for extremely small and unphysical pertur- 
bations (or times). For instance, according to the argument by van Kampen, 
if one wants that the FDR holds up to Is when applied to the electrons in a 
conductor then a perturbing electric field should be smaller than 10~^'^V/m, in 
clear disagreement with the experience. 

This result, at first glance, sounds as bad news for the statistical mechanics. 
However this criticism had the merit to stimulate for a deeper understanding 
of the basic ingreedients for the validity of the FDR relation and its validity 
range. On the other hand the success of the linear theory, for the computation 
of the transport coefficients (e.g. electric conductivity) in terms of correlation 
function of the unperturbed systems, is transparent, and its validity has been, 
directly and undirectly, tested in a huge number of cases. Therefore it is really 
difficult to believe that the FDR cannot be applied for physically relevant value 
of the perturbations. 

Kubo suggested that the origin of this effectiveness of the Linear-Response 
theory may reside in the "constructive role of chaos" because, as Kubo sug- 
gested, ^^instability [of the trajectories] instead favors the stability of distribution 
functions, working as the cause of the mixing" |20| . The following derivation [17] 
of a generalized FDR supports this intuition. 

2.2 Generalized FDR for stationary systems 

Let us briefly discuss a derivation, under rather general hypothesis, of a gen- 
eralized FDR (also in non equilibrium or non Hamiltonian systems), for which 
the van Kampen critique does not hold. Consider a dynamical system x(0) — >■ 
x(t) = t/*x(0) with states x belonging to a A^-dimensional vector space. For 
the sake of generality, we will consider the case in which the time evolution 
can also be not completely deterministic {e.g. stochastic differential equations). 
We assume the existence of an invariant probability distribution jo(x), for which 



some "absolutely continuity" type conditions are required (see later), and the 
mixing character of the system (from which its ergodicity follows) . Note that 
no assumption is made on N. 

Our aim is to express the average response of a generic observable A to a 
perturbation, in terms of suitable correlation functions, computed according to 
the invariant measure of the unperturbed system. At the first step we study 
the behaviour of one component of x, say a;^, when the system, described by 
p(x), is subjected to an initial (non-random) perturbation such that x(0) — >■ 
x(0) + Axq q This instantaneous kick modifies the density of the system into 
yo'(x), related to the invariant distribution by /c'(x) = p(x — Axq). We introduce 
the probability of transition from Xq at time to x at time t, w(xo,0 — >■ x, i). 
For a deterministic system, with evolution law x(i) = f7*x(0), the probability 
of transition reduces to ^(xo, — >■ x, i) = S{x — C/*Xo), where S{-) is the Dirac's 
delta. Then we can write an expression for the mean value of the variable a;^, 
computed with the density of the perturbed system: 

(xi{t)\ = x.,p'{xo)w{xo,0^x,t)dxdxo . (8) 

The mean value of Xi during the unperturbed evolution can be written in a 
similar way: 

(x^{t)\ = Xjp(xo)u;(xo,0 ^ x,t)(ixdxo . (9) 

Therefore, defining Sxi = (xi)' — (xi), we have: 

Sxi{t) = / Xi — p(xo)w(xo,0 ^^ X, i)dxdxo 

= (a-,(t)F(xo,Axo)) (10) 

where 



^(xo,Axo) = 



p(xo - Axo) - /9(xo) 



(11) 



p(xo) 

Let us note here that the mixing property of the system is required so that the 
decay to zero of the time-correlation functions assures the switching off of the 
deviations from equilibrium. 

For an infinitesimal perturbation Sx{0) = {Sxi(0) ■ ■ ■ <5a;jv(0)), if p(x) is non- 
vanishing and differentiable, the function in ([TT|) can be expanded to first order 
and one obtains: 



Sx^{t) = -^(a;,(t) 



dlnp{x) 



dxj 



)Sx,{0) 



t=Oj 



^i?,,,(t)fa,(0) (12) 



^The study of an "impulsive" perturbation is not a limitation, e.g. in the linear regime 
from the (differential) linear response one understands the effect of a generic perturbation. 



which defines the linear response 

d\np{x) 



i?,j(t) = - x,(t) 



''''^' '^''^' dx, 



(13) 

t=o/ 



of the variable Xi with respect to a perturbation of Xj . 

We note that in the above derivation of the FDR relation we never used any 
approximation on the evolution of (5x(i). Starting with the exact expression ([9]) 
for the response, only a linearization on the initial time perturbed density is 
needed, and this implies nothing but the smallness of the initial perturbation. 
We have to stress again that, from the evolution of the trajectories difference, 
one can define the leading Lyapunov exponent A by considering the absolute 
values of (5x(i): at small |(5x(0)| and large enough t one has 

/in |5x(t)|\ ~ In |,5x(0)| + \t . (14) 

On the other hand, in the FDR issue one deals with averages of quantities with 
sign, such as (Sx(i). This apparently marginal difference is very important and 
it is at the basis of the possibility to derive the FDR relation avoiding the van 
Kampen's objection. 

2.3 Remarks on the invariant measure 

At this point one could object that in a chaotic deterministic dissipative system 
the above machinery cannot be applied, because the invariant measure is not 
smooth at alo In chaotic dissipative systems the invariant measure is singular, 
however the previous derivation of the FDR relation is still valid if one considers 
perturbations along the expanding directions. A general response function has 
two contributions, corresponding respectively to the expanding (unstable) and 
the contracting (stable) directions of the dynamics. The first contribution can 
be associated to some correlation function of the dynamics on the attractor 
(i.e. the unperturbed system). On the contrary this is not true for the second 
contribution (from the contracting directions), this part to the response is very 
difficult to extract numerically [21]. Nevertheless a small amount of noise, that 
is always present in a physical system, smoothes the p(x) and the FDR relation 
can be derived. We recall that this "beneficial" noise has the important role of 
selecting the natural measure, and, in the numerical experiments, it is provided 
by the round-off errors of the computer. We stress that the assumption on 
the smoothness of the invariant measure allows one to avoid subtle technical 
difficulties. 

In Hamiltonian systems, taking the canonical ensemble as the equilibrium 
distribution, one has Inp = — /37J(Q,P) + constant. Recalling Eq. (|T3)) . if we 
indicate there by Xi the component q^ of the position vector Q and by Xj 
the corresponding component pk of the momentum P, from the Hamilton's 



■^Typically the invariant measure of a chaotic attractor has a multifractal character and its 
Rcnyi dimensions dq are not constant. 



equations {dqk/dt ~ dH/dpk) one has the differential form of the usual FDR [H 
[20] 

^k{t) J ,,,dqk{0)\ d/ 



Spk{0) 



/3(g.(i)^)=-/3|(%(t)%(0)). (15) 



When the stationary state is described by a multivariate Gaussian distribu- 
tion, 

Inp(x) = — - y^ctijXiXj + const. (16) 



2 



with {aij} a positive matrix, and the elements of the linear response matrix can 
be written in terms of the usual correlation functions: 



^«.j(*) = 5I"^^fe(^'W^fe('^)) ■ (1^) 



The above result is nothing but the Onsager regression originally obtained for 
linear Langevin equations. It is interesting to note that there are important non- 
linear systems with a Gaussian invariant distribution, e.g. the inviscid hydro- 
dynamics |22U23| . where the Liouville theorem holds, and a quadratic invariant 
exists. 

In non Hamiltonian systems, where usually the shape of p(x) is not known, 
relation (|13p does not give a very detailed quantitative information. However 
from it one can see the existence of a connection between the mean response 
function Rij and a suitable correlation function, computed in the non perturbed 
systems: 

(x,(f)/,(x(0))) , with /^.(x) = -^, (18) 

where, in the general case, the function fj is unknown. 

Let us stress that in spite of the technical difficulty for the determination of 
the function fj, which depends on the invariant measure, a FDR always holds 
in mixing systems whose invariant measure is "smooth enough" . We note that 
the nature of the statistical steady state (either equilibrium, or non equilibrium) 
has no role at all for the validity of the FDR. 

We close this Section noting that, as well clear even in the case of Gaussian 
variables, the knoweledge of a marginal distribution 

Pi{xi)= p{xi,X2,--.-)Y\_d'Xj (19) 

is not enough for the computation of the auto response: 



RiAW r -( ^iv) 



dxi 



) , (20) 

t=o/ 



This observation is particularly relevant for what follows. Consider that the de- 
scription of our system has been restricted to a single slow degree of freedom, for 



instance the coordinate of a colloidal particle ruled by the Langevin equation. 
The above discussion has indeed made clear that there can be other fluctuating 
variables coupled to the one we are interested in, and it is not correct to project 
them out by the marginalization in Eq. (J19p . Conversely, a stationary propa- 
bility distribution with new variables coupled to the colloidal particle position 
must be taken into account. In the examples discussed in the last two sections, [5] 
and [6] these additional fluctuating variables turn out to be local fields coupled 
with the probe particle. 

2.4 Generalized FDR for Markovian systems 

In the previous section we have discussed a formula useful to write an explicit 
expression for the response to an external perturbation at stationarity in a case 
where the invariant measure is known. The possibility to work out explicitly a 
response formula is indeed quite general. In particular, for Markov processes a 
general FDR has been derived [24l [25l [26] , and also experimentally verified p7l 
[28], which also holds for non-stationary, aging processes, even in the absence of 
detailed balance. Here we report a derivation which strictly follows the one given 
in |24| . An interesting outcome for the purpose of the whole discussion will be 
the evidence that the response function includes quite generally the correlation 
of the field of interest with a local field. 

Let us briefly recall that a Markov process is univocally defined by the initial 
distribution probability p(x, 0) and the transition rates W{'x! — >■ x) from state 
x' to state x, with normalization 

Y^ VK(x' ^ x) = 0. (21) 

X 

The two-time conditional probability is then approximated by 

P(x, t + A<|x', t) = 5x.x' + VK(x' -^ x)At + 0(Af2), (22) 

where At is a small time interval. The average of a given observable A(x) at 
time t reads 

{A{t)) = J2 ^(x)P(x, t|x', t')p(x', t'). (23) 

x.x' 

Now we discuss a perturbation to our system in the form of a time-dependent 
external field h{s), which couples to the potential V{x) and changes the energy of 
the system from ?^(x) to ?^(x) — h{s)V{x). The dynamics in the presence of the 
pertubation is decribed by a new Markov process defined through the perturbed 
transition rates W^''(x|x'). There is not a uni vocal prescription constraining on 
the choice of a particular form of W'^. Here we focus on the FDR following 
from the particular choice which obeys the so-called local detailed balance |29| . 
In this case, to linear order in h, the perturbed transition rates read, for x 7^ x'. 



W''{x' ^ x) = W{x' ^ x) |l ~ ^ [y(x) - \/(x')] + Af(x,x')| , 



(24) 



where (3 is the inverse temperature and M(x, x') is an arbitrary function of 
order /3h, symmetric with respect the exchange of its arguments. The diagonal 
elements are obtained imposing the normalization condition (|2ip . Notice that, 
once the local detailed balance is imposed, there remains a further degree of 
arbitrariness in the choice of the function M. The dependence of the FDR on 
the form of such function is studied in [211 HU • Here we focus on the particular 
case M = 0. 

Let us consider an impulsive perturbation turned on at time s for the time 
interval At. The linear response function of the observable A is, a.t t > s + At, 
with t - s > At, 

1 6{A{t))h 



R{t,s) 



lim 



(25) 



h=0 



At-i-o At Sh{s) 

where (. . .)h denotes an average on the perturbed dynamics. The derivative 
with respect to h is written in terms of the conditional probabilities as follows 



5{Ait))H 



Shis) 



h=0 



E 

x,x' x" 



A(x)P(x,t|x',.s+At) 



(5P''(x',.s + At|x",.s) 



Sh 



P(x",s) 



h=0 



(26) 
In order to derive a FDR relating the response function to correlation func- 
tions computed in the unperturbed dynamics, the derivative with respect to the 
external field in the rhs of Eq. ([26]) has to be worked out explicitly. To do this, 
we notice that, enforcing the normalization condition (|2ip . the transition rates 
can be separated in diagonal and off-diagonal contributions, namely 



M/"(x' ^ x) = -<5x,x' Y. W^"(x'^x") + (l-<5x,x')W^"(x'^x). (27) 

x"^x' 

Then, using Eqs. ([22]), ^ and ^, one obtains 



1 SP^{x,s + At\yi',s) 



At 



5h 



h=0 



f ^^x,x' J2 M^(x"|x')[X^(x')-y(x")] + (l-5x,x')W^(x'^x)[F(x)-F(x')] 

x"7^x' 



(28) 



Such separation allows us to single out two different terms. Indeed, substituting 
Eq. ^ into Eq. ^, one obtains 
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S{Ait)), 



Sh{s) 



P 



h=0 



J2 ^WP(x, i|x', s + At)(AtW{x' -^ x")[V^(x') - F(x")]P(x', s) 



+ AiVK(x" ^ x')[F(x') - T/(x")]P(x",s)}. 



Then, exploiting the time translation invariancc of the conditional probability, 
namely P(x, i + Ai|x',i) = P(x, t|x',i — At), in the first line of the above 
equation, and using Eq. ((22|) in the second one, in the limit Ai — >■ one obtains 
the response function 



R(t,s) 



/3 



d{A{t)V{s)) 
ds 



{A(t)B{s)) 



where 



Bis) EE B[^is)] = ^{y(x") - V[^{s)]}WMs) ^ x"] 



(30) 
(31) 



is an observable quantity, namely depends only on the state of the system at 
a given time. The relation ()30p is known since a long time in the context of 
overdampcd Langcvin equation for continuous variables |30| . 

For instance, in the case of the overdampcd Langevin dynamics of a colloidal 
particle diffusing in a space dependent potential U{x), 



x{t) 



dUjx) 
dx 



2TC(t), 



(32) 



with C,{t) a white noise with zero mean and unit variance, one can derive a 
response formula, with respect to a perturbing force F, analogous to Eq. (|30|) 
which reads as: 



S{<t))i 
6F{s) 



d{x{t)x{s)) 



ds 



{xit)B[x{s)]) 



(33) 



with B[x{s)] ~ —dU / dx\x(s)- At equilibrium it can be easily proved that 
{x{t)B[x{s)]) = ~d{x{t)x{s)) /ds, recovering the standard FDT formula. Dif- 
ferently, out of equilibrium the contribution coming from the local field B[x{s)] 
must be explicitly taken into account, as will be shown in examples in the next 
sections. 



(29) 



3 Random walk on a comb lattice 

3.1 Anomalous diffusion and FDR 

As discussed above for Brownian motion, in the absence of external forcing one 
has, for large times i — > oo. 



(x(i))=0 , {x\t))^2Dt, 



(34) 
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where x is the position of the Brownian particle and the average is taken over 
the unperturbed dynamic. Once a smah constant external force F is applied 
one has a linear drift 

fe(i) = {x{t))F - {x{t)) ~ ^iFt (35) 

where {■ ■ ■)f indicates the average on the perturbed system, and fi is the mobility 
of the colloidal particle. It is remarkable that (a;^(i)) is proportional to Sx{t) at 
any time: 

a^ = A, (36) 

Sx{t) PF' ^ ' 

and the Einstein relation (see Eq. S]) holds. 

On the other hand it is now well established that beyond the standard dif- 
fusion, as in (j34p . one can have systems with anomalous diffusion (sec for in- 
stance [311 lia EH [Ml EH [36]), ie. 

{x'^{t))^f' with iy^l/2. (37) 

Formally this corresponds to have D = oo ii ly > 1/2 (supcrdiffusion) and D = Q 
\i V < 1/2 (subdifFusion). In the following we will limit the study to the case 
V < 1/2. 

It is quite natural to wonder if (and how) the FDR changes in the pres- 
ence of anomalous diffusion, i.e. if instead of (|34|) . Eq. ([37]) holds. In some 
systems it has been shown that ([36]) holds even in the subdiffusive case. This 
has been explicitly proved in systems described by a fractional-Fokker-Planck 
equation [S^, see also [38l [39]. In addition there is clear analytical [40] and 
numerical [41] evidences that ([36]) is valid for the elastic single file, i.e. a gas 
of hard rods on a ring with elastic collisions, driven by an external thermostat, 
which exhibits subdiffusive behavior, (x^) ~ i^/^ |42] . 

Here we discuss the validity of the fluctuation-dissipation relation in the 
form pop for a system with anomalous diffusion which is not fully described 
by a fractional Fokkcr-Planck equation. In particular we will investigate the 
relevance of the anomalous diffusion, the presence of non equilibrium conditions 
and the (possible) role of finite size. The example discussed here is the diffusion 
of a particle on a comb lattice. The dynamics of the model is defined by tran- 
sition rates and therefore a straightforward application of the generalized FDR 
introduced in Sections |2.4[ Eq ([30]) . is possible. 

3.2 The transition rates of the model 

The comb lattice is a discrete structure consisting of an infinite linear chain 
(backbone) , the sites of which are connected with other linear chains (teeth) of 
length L [43]. We denote by a; G (— oo, oo) the position of the particle performing 
the random walk along the backbone and with y € [—L,L] that along a tooth. 
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Figure 1: Left panel: {x'^{t))o/L vs t/L"^ is plotted for several values of L in the 
comb model. Right panel: {x'^{t))o and the response function Sx{t) for L = 512. 
In the inset the parametric plot dx{t) vs {x^{t))o is shown. 



The transition probabilities from state x ^ {x, y) to x' = {x' , y') are: 



W'^[{x,0)^ (a;±l,0)] 

W^'^[(x,0)-^(x,±l)] 

W''[{x,y)^{x,y±l)] 



1/4 ±d 

1/4 

1/2 ioTy^O,±L. 



(38) 



On the boundaries of each tooth, y = ±L, the particle is reflected with proba- 
bility 1. The case L = oo is obtained in numerical simulations by letting the y 
coordinate increase without boundaries. Here we consider a discrete time pro- 
cess and, of course, the normalization J2(x' y') ^'^[(^^y) ~^ {^' iV')] = 1 holds. 
The parameter d S [0,1/4] allows us to consider also the case where a constant 
external field is applied along the x axis, producing a non zero drift of the par- 
ticle. A state with a non zero drift can be considered as a perturbed state (in 
that case we denote the perturbing field by e), or it can be itself the stationary 
state where a further perturbation can be added changing d — > d + e. 



3.3 Anomalous dynamics 

Let us start by considering as a reference state the case d = 0. For finite teeth 
length L < oo, we have numerical evidence of a dynamical crossover, at a time 
t*, from a subdiffusive to a simple diffusive asymptotic behaviour (see Fig. [ij 



'm 



Ct^/^ t<t*{L) 
2D{L)t t>t*{L), 



where C is a constant and D{L) is an effective diffusion coefficient depending 
on L. The symbol (. . .)o denotes an average over different realizations of the 
dynamics ((38)) with d = and initial condition a;(0) = y(0) = 0. We find 
t*{L) ~ L^ and D{L) - 1/L. In the left panel of Fig. [1] we plot {x'^{t))o/L as 
function of t/i^ for several values of L, showing an excellent data collapse. 
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In the limit of infinite teeth, L — >■ cxj, D — > and t* — > oo and the system 
shows a pure subdiffusive behaviour j44j 



(X 



2(i))o^ti/2. (39) 



In this case, the probability distribution function behaves as 

( Ni V^ 
Po{x,t)^t-^'^e~'y7T^) , (40) 

where c is a constant, in agreement with an argument a la Flory [31] . The 
behaviour (PH)) also holds in the case of finite L, provided that t < t* . For larger 
times a simple Gaussian distribution is observed. 

In the comb model with infinite teeth, the FDR in its standard form is 
fulfilled, namely if we apply a constant perturbation e pulling the particles 
along the f-d lattice one has numerical evidence that 

(z2(t))o~CMt)~ti/2^ (41) 

where 6x is the average in the presence of the perturbation e. Moreover, the 
proportionality between {x'^{t))o and Sx{t) is fulfilled also with L < oo, where 
both the mean square displacement (MSD) and the drift with an applied force 
exhibit the same crossover from subdiffusive, ^ t^'^, to diffusive behavior, ~ t 
(see Fig.[Tl right panel). Therefore what we can say is that the FDR is somehow 
"blind" to the dynamical crossover experienced by the system. When the per- 
turbation is applied to a state without any current, the proportionality between 
response and correlation holds despite anomalous transport phenomena. 

3.4 Application of the generalized FDR 

Our aim here is to show that, differently from what depicted above about the 
zero current situation, within a state with a non zero drift the emergence of 
a dynamical crossover is connected to the breaking of the Einstein FDR (|3]). 
Indeed, the MSD in the presence of a non zero current, even with L = oo, shows 
a dynamical crossover 

{x'it))d^ae^'' + bt, (42) 

{[xit)-{x{t))a]')dr^a't'/^ + b't, (43) 

where a, b, a' and b' are constants, whereas 

I^4t)-t'^\ (44) 

with Sxd{t) = {x{t))d+e — {x{t))d'- hence at large times the Einstein rela- 
tion breaks down (see Fig. [5]). Notice that the proportionality between re- 
sponse and fluctuations cannot be recovered by simply replacing {x'^(t))d with 
{x^{t))d — {x{t))'^, as it happens for Gaussian processes (see discussion below). 
The constants a' and b' can be computed analitically in the case L = oo [45) . 
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The first moment {x{t))d+e of Pd+e{x,t) is not proportional to the second 
cumulant of Pd{x, t), namely {x^{t))d — {x{t))^. In order to point out a relation 
between such quantities, we need a generalized fluctuation-dissipation relation. 

According to the definition (|38p . one has for the backbone 

(46) 

where W^ = 1/4, and the last expression holds under the condition d/W'^ ^ 1. 
The above expression can be rewritten in the form of a local detailed balance, 
Eq. ([24]), with F(x) = x and /3 = 1, yielding, for {x,y) ^ {x\y'), 

W+^x^y) ^ {x\y')] = Wlix^y) ^ {x\y')]e'^'^^'-^\ (47) 

where /i(e) = 2e/M^°. Using Eqs. ^^ and ([31]) we obtain 

50d_ ^ (^W)-^+^^^- (^W)'^ = 1 [(O(t)x(t)). - (O(t)x(O)). - (O(t)A(t,0)).] , 

(48) 
where O is a generic observable, and A{t, 0) = X]t'=o ^(^')i with, in this case, 

B[{x,y)] = ^ (a;'-x)M/^[(x,y) -^ {x\y')]. (49) 

(2;',j/') 

Recalling the definitions (|38| . from the above equation we have B[(a;,y)] = 
2dSyfi and therefore the sum on B has an intuitive meaning: it counts the time 
spent by the particle on the x axis. The results described in section [231 can be 
then read in the light of the fluctuation-dissipation relation (|48l) : 

i) Putting 0{t) = x{t), in the case without drift, i.e. d = 0, one has B = 
and, recalling the choice of the initial condition a;(0) = 0, 

fe {x{t)), - {x{t))o _ 1 , 2 



Me) his) -2<^(^)>- (^°) 

This explains the observed behaviour (|4ip even in the anomalous regime and 
predicts the correct proportionality factor, 6x{t) = e/W'^{x'^{t))o. 
ii) Putting 0{t) = x{t), in the case with d 7^ 0, one has 



h{e) 2 



[{x'it))d^{x{t)A{t,0))d]. (51) 



This explains the observed behaviours (|^^ and ([33]): the leading behavior 
at large times of {x^{t))d ~ t, turns out to be exactly canceled by the term 
(x(i)A(i, 0))d, so that the relation between response and unperturbed correla- 
tion functions is recovered (see Fig. [2]). 
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Figure 2: Response function (black line), MSD (red dotted line) and second 
cumulant (black dotted line) measured in the the comb model with L = oo, field 
d = 0.01 and perturbation e — 0.002. The correlation involving the quantity 
B (green dotted line) yields the right correction to recover the full response 
function (blue dotted line), in agreement with the FDR 



iii) As discussed above, it is not enough to substitute {x'^{t))d with {x'^{t))ii — 
(a;(i))^ to recover the proportionality with Sxd{t) when the process is not Gaus- 
sian. This can be explained in the following manner. By making use of the 
second order out of equilibrium FDR derived by Lippiello et al. in [46l |47l |48] , 
which is needed due to the vanishing of the first order term for symmetry, we 
can explicitly evaluate 



{x\t))a = {x\t)), + h\d)]- 



i(x^(t))o + i(x^(t)A(2'(i,0))c 



(52) 



where ^(2) (i, 0) = Y.\'=o B^^^ (f) with ^(2) = - J^x' ix'-x)^W[{x, y) -> (x', y')] 
— l/25yfl. Then, recalling Eq. ((50)) . we obtain 



{x\t))d-{x{t))l 
{xHt))o + h^{d) 



l{a:\t))o + l{xHt)A^'\t,0))o-\{xHt))l 



(53) 



Numerical simulations show that the term in the square brackets grows like t 
yielding a scaling behaviour with time consistent with Eq. (j43p . On the other 
hand, in the case of the simple random walk, one has i?^^) = — i and A^'^'>(t, 0) = 



16 



-t and then 



X^t))o + h^d) 



l{x^{t))o-lt{x'{t))o-l{x'it))l 



(54) 



Since in the Gaussian case (x^(t))o — 3(x^(i))o and (x^(t))o — t, the term in the 
square brackets vanishes identically and that explains why, in the presence of a 
drift, the second cumulant grows exactly as the second moment with no drift. 

4 Entropy Production 

Non equilibrium regimes are always characterized by some sort of current flowing 
across the systems. In the following sections we will suggest, showing examples, 
some connections between such non-equilibrium currents and the new degrees 
of freedom entering the fluctuation dissipation relation, while in the present one 
we introduce a general symmetry relation of the probability distribution of such 
currents, the Fluctuation relation. 

In general, when one deals with non-equilibrium dynamics, very few results 
independent from the details of the model are available. Actually in the last 
decade a group of relations, known with the name of "fluctuation relations" have 
captured the interest of the scientific community, especially for the vast range of 
applicability. Initially, a numerical evidence given by Evans and collaborators 
showed a particular symmetry in the distribution of an observable of a molecular 
fluid under shear. In a second moment, such a symmetry has been proved 
as a theorem by Gallavotti and Cohen, under quite general hypothesis. The 
interested reader can see, among others [33] ■ For a review of large deviation 
theory applied to non-equilibrium systems, see [50] . 

According to the point of view here adopted we will focus on systems in which 
some randomness is present. Thanks to this assumption, it is possible to skip 
several technical problems and some forms of fluctuation theorems for stochastic 
systems can be used. Among others, we are going to use the Lebowitz-Spohn 
functional, valid for Markovian dynamics. In order to fix ideas, let us consider 
a one-dimensional process discrete in time (the generalization to the other cases 
is straightforward) and let us identify a trajectory in a time interval [0,t]: 

a = {x(0),x(l),...,a:(i)} (55) 

Clearly, because of the stochastic nature of the dynamics it is possible to asso- 
ciate to the trajectories a probability P{Ut); using the Markovian nature of the 
process one has 

t 
Piflt) = p{x{0}) Y[ W{x{?i - 1) ^ x{n)) (56) 

n=l 

Analogously, one can consider the time-reversed trajectory, i.e. fit = x{t), x{t ~ 1), . . . , x(0) 
with its probability P{flt). For clarity, here we are considering variables which 
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are even in time. This is not the most general case: for example, when one 
deals also with velocities one must take into account the different parity of the 
variables. Then the Lebowitz-Spohn functional is defined as 

Note that this quantity, in general, is not related to a specific thermodynamic 
observable. However, in what follows, we will call it "entropy production" ac- 
cording to the recent literature. If the dynamics satisfies detailed balance con- 
ditions, one has that P{flt) = P{^t) and then St is identically equal to zero. 
On the contrary, when the detailed balance condition is broken its average is 
strictly positive and does increase in time 

(Et) - i for i large (58) 

where the brackets mean an average on the space of trajectories in the stationary 
ensemble. In this sense, one of the main features of this quantity is that it 
captures the "non-equilibrium" nature of the system. 

Let us discuss with a pedagogical example how the entropy production is 
related to non-equilibrium currents. Consider a Markov process where the per- 
turbation of an external force F inducing non-equilibrium currents enters in the 
transition rates according to local detailed balance condition 

Wf{x'-^x) Wo{x'^x) ' ^ ' 

with j{x — ^ x') the current associated to the transition x — > x' , which obey the 
symmetry property j{x — >■ x') = —j{x' — >■ x). According to the definition of 
entropy production ([57)) one finds, for large times, 

— = 2I3F- Y^ J{x{n - 1) ^ x(n)) = 2(3FJ{t), (60) 

where J{t) is the average current over a time window of duration t. The fluc- 
tuation theorem is a symmetry property of the probability distribution of the 
variable y — Y,t/t which reads: 



P{y) .._ P(2/3Fj(i)) 



P{-y) Fi-2(3FJ{t)) 



e 



2l3FJ{t] 



(61) 



Namely the fluctuation theorem describes a symmetry in the fluctuations of 
currents. Also, for large times we can assume a large deviation hypothesis 
P{y) ^ e~*^^y\ with S{y) a Cramer function. For small fluctuations around the 
mean value oiy the Cramer function can be approximated to S{y) = S{2pFJ) ~ 
I3^F^{J — JY /a'^j, where aj is the variance. The fluctuation relation reads as 
S{y) — S{—y) ~ y; in the Gaussian limit {y close to y) the previous constraint 
can be easily demonstrated to be equivalent to J /F = I3aj, which is nothing 
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but the standard fluctuation dissipation relation. Therefore the fluctuation re- 
lation, which in some simple situation can be directly related to the fluctuation- 
dissipation relation, is a more general symmetry to which we expect to obey 
the fluctuating entropy production. For a more general discussion of the link 
between the Lebowitz-Spohn entropy production and currents, see [511 [S^- The 
remarkable fact appearing in equation (j6ip is that it does not contain any free 
parameter, and so, in this sense, is model- independent. 

5 Langevin processes without detailed balance 

Sometimes the properties of statistical systems are investigated studying the 
diffusional properties of probe particles which customarily have larger mass and 
size compared with the constituents of the environment. This is for instance 
the case of Brownian motion, in which the erratic motion of flower-dust within 
water molecules was considered. The central point is that such probe particle 
is always coupled to a small portion of the whole system, so that small scale 
fluctuations always do matter. We are going to study the feedback of local 
fluctuations of the fluid surrounding the probe on the dynamics of the probe 
itself when the environment is out of equilibrium. The dynamics of a single 
probe particle in contact with an equilibrium thermal bath is well described by 
a single stochastic differential equation: 

MV{t) = - / dt' T{t - t')V{t') + £{t'), (62) 



where £{t) is a gaussian white noise with zero mean and variance {£{t)£(t')) = 
r(|t — t'l) (second kind FDT). The equilibrium condition is always guaranteed 
by the proportionality between memory kernel and noise autocorrelation. The 
simpler way to account for the coupling between our probe particle with non 
equilibrium fluctuations of the environments is by breaking the second kind 
FDT. 

In the rest of the section, we will discuss the following points: 

• how non-equilibrium arises in a multidimensional linear Langevin model; 

• how the coupling between the variable of interest and others do matter in 
the response formula; 

• what form is taken by the Fluctuation-Relation in such a system. 

5.1 Markovian linear system 

The coupling between the probe and the portion of fluid surrounding it, Eq. (|62p . 
in certain regimes may be effectively described in terms of two coupled Langevin 
equations: one variable is the velocity of our tagged particle and the other 
is a local force field. Such a simplification can be realized in not too dense 
cases, where the memory kernel F has a single characteristic finite time: a more 
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detailed discussion can be found in [S31[S11 and a specific example, in the context 
of granualr systems, is given in section |6l The system can be then put in the 
following form |55) : 

MV = -T{V -U) + V2rTi(^i (63) 

M'ii = -T'u ~rv + y/2rT2(i32, 

where M and M' are "masses" , T and F' arc sort of viscosities and Ti and T2 
are two different "energy scales" . 

Model (|63p . in a more compact form, reads 

f.-AX + ^, (64) 

where X = (Xi,X2) e (p = {(l>i,(l>2) are 2-dimensional vectors and A is a real 
2x2 matrix and (j)(t) a Gaussian process, with zero mean and covariance matrix: 

{Mt')hit))^2Sit~t')D,^, (65) 

and 

A=( t -^¥) D=\ ^;: r'T. ). (66) 




The stability conditions on the dynamical matrix which guarantee the existence 
of a stationary state are Tr(A) > and det(^) > which are clearly fulfilled in 
the present case. The invariant measure of the steady state is represented by a 
2-variate Gaussian distribution: 

p{X)=Ncxp(-^Xa-'x\ (67) 

where A^ is a normalization coefficient and the matrix of covarianccs a is ob- 
tained by solving 

D = Aa + aA^, (68) 

which yields: 



^ + eAT GAT 



M 



GAT ^ + feAT 



(69) 



where 8 = (Y+T')(M'r+Mr') ^^"^ ^-^ ~ Ti — T2. It is now well clear that 
when Ti = T2 the two variables are uncorrelated. This is the fingerprint of an 
equilibrium condition, as we shall see in the following sections. 

5.2 Fluctuation-response relation 

We can now explicitly study the fluctuation and response properties of the 
system since the dynamics is linear. First, the correlation matrix Cij(t,t') = 
{Xi{t)Xj{t')), in the stationary state is time-translational invariant, i.e. it only 
depends on the time difference t — t' . Then, using the equation of motion, it is 
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immediate to verify that C{t) = —AC{t), with initial condition given by the co- 
variance matrix, i.e. C(0) = a. The corresponding solution (in the matrix form) 
is: C{t) — e~^*a Note that, in general a and A do not commute. Moreover, 
considering the gaussian shape of the steady state distribution function (p7)) it 
is possible to calculate the response function of the system, by a straightforward 
application of ([T7|) : 



^^^^ " WW) = E^y'(^^W^^(0)) = C(t)a-i. (70) 

In the case of our interest, by considering a perturbation applied to the variable 
V, one obtains: 



^=a^,\V{t)Vm+a^,\V{t)Um. (71) 

As it appears clear from equation ([7T]) , the response of the variable V to an ex- 
ternal perturbation in general is not proportional to the unperturbed autocorre- 
lation. It is possible to observe the following scenario: in the case Ti = T2 = T, 
i.e. AT = 0, (7 is diagonal with an = -^ and 0-22 = -jp, independently of the 
values of other parameters, a direct proportionality between Cvv and Rvv is 
obtained. This is not the only case for this to happen: also for F ^ F', or 
F ^ F', (712 goes to zero and a direct proportionality between response and 
autocorrelation of velocity is recovered. More in general, when T2 ^ Ti, the 
coupling term (T12 differs from zero and a "violation" of the Einstein relation 
emerges, triggered by the coupling between different degrees of freedom. 

The situation is clarified by Fig. [3j where the Einstein Relation is violated 
and the generalized FDR holds: response Rvv{t), when plotted against Cyy(i), 
shows a non-linear (and non-monotonic) graph. Anyway a simple linear plot is 
restored when the response is plotted against the linear combination of correla- 
tions indicated by formula ([7T|) . In this case it is evident that the "violation" 
cannot be interpreted by means of any effective temperature |56j : on the con- 
trary it is a consequence of having "missed" the coupling between variables V 
and U, which gives an additive contribution to the response of V. 

The above consideration can be easily generalized to the case where several 
auxiliary fields are present: let us suppose that there are N — 1 variables at 
temperature Ti (the velocity V and N — 2 auxiliary variables Ui) and one at 
temperature T2. In this case it is possible to show that the off-diagonal terms 
in a arc proportional to (Ti — T2) and a qualitative similar analysis can be 
performed. 

It is useful to stress the role of Markovianity, which is relevant for a correct 
prediction of the response. In fact the marginal probability distribution of 
velocity Pm{V) can be computed straightforwardly and has always a Gaussian 
shape. By that, one could be tempted to conclude, inserting Pm{V) inside , that 
proportionality between response and correlation holds also if T2 7^ Ti. This 
conclusion is wrong, as shown in equation ((20| because the process is Markovian 
only if both the variable V and the "hidden" variable U are considered. 
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Figure 3: Free particle with viscosity and memory, whose dynamics is given 
by Eq. (j63|) : here we show the parametric plot of the velocity response to an 
impulsive perturbation at time 0, versus two different correlations. The Einstein 
relation, which is not satisfied, would correspond to a linear shape with slope 
1 for the red dashed line. The black line shows that the Generalized response 
relation holds. 



Wc conclude by mentioning a sentence of Onsager: how do you know you 
have taken enough variables, for it to be Markovian? |57] . We have seen in this 
section how a correct Response-Fluctuation analisys is intimately related to this 
important caveat. 

5.3 Entropy production 

Due to the linearity of the model ((64|). the explicit calculation of the entropy 
production according to the Lebowitz-Spohn definition is quite straightforward. 
There is a one to one mapping between the trajectories of the system and the 
realizations of the noise which allows to write the probability of each trajectory 
as a path-integral over all the possible realizations of noise: 



P{{Ms)Yo) = / Vcj,P{cf,)d{± + AX - , 



(72) 



Due to the gaussianity of P(0) the above integral is easily evaluated yielding 
the well know Onsager-Machlup expression for the path probabilities: 



P({X(s)}*) ^ e-(^+^^)^"(^+^^). 



(73) 



In our model it is possible to compute the entropy production functional of 
a single trajectory. Let us consider a general trajectory {X(s)}g and its time- 
reversal {X(s)}g. Lebowitz and Spohn defined a fiuctuating entropy production 
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functional Wt as follows j52j : 

p[X(0)]P({X(,s)}^) 
* %[X(t)]P({X(.)}*) ^ ^ 

with 

bt = log{/[X(0)]} - log{/[(X)]}, (75) 

where p[X(0)] is the stationary distribution, i.e. a bivariate gaussian with co- 
variance given by equation (j68p . Lebowitz and Sphon have shown that the 
average (over the steady ensemble) of Wt, if detailed balance is not satisfied, 
increases linearly with time, while the term bt, usually known as "border term", 
is usually negligible for large times, unless particular conditions of "singularity" 
occur [SillMllSO]- 

For semplicity of the notation, let us define ( t^ /v\ ) — ^ ^X. In order 

V -f'2(X) J 

to write down an explicit expression for the entropy production, it is necessary 

to establish parity of variables under time-reversal (i.e. positions are even and 

velocities are odd under time inversion transformation). Let us assume that 

under time reversal Xi = eiXi, where e^ can be +1 or —1, using also eX = 

(€1X1,62X2). Then it is possible to define 

Frm = Im^) - e,F,{eX)] ^ -e,FrieX) (76) 

Fr(X) = ^[F,(X} + e,F,{eX)]^e,Fr{eX). (77) 

Given this notation |61| it is possible to write down a compact form for the en- 
tropy productiorO by simply substituting equation ([73| into (|74|) and obtaining: 



S* - E ^fefe / d'P^' \^^ - ^^^'1 ' (78) 



where the sum is over k such that D^k 7^ 0. Formula (|78p is valid also in presence 
on non-linear terms and with several variables |54| . 

For the model system in Eq. (|B5|. it appears clear that the variable V, being 
velocity, is odd under the change of time, while variable C/, which represents 
the local force field (divided by F), is even. According to the above formula the 
entropy production reads as: 

^tc^r(^^^-l-^ jyit')uit')dt'. (79) 

Note that, from (|79p . one recovers the equilibrium case, i.e. Ti = T2, for 
which there is no entropy production. Remarkably, the equality of the two 
temperatures is the same equilibrium condition given by the FDT analysis. 



•'Note that if the correlation matrix of the noise is not diagonal, this formula is slightly 
different. 
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6 Granular intruder 

Models of granular fluids |62| are an interesting framework where the issues 
discussed in the previous sections can be addressed. Due to dissipative interac- 
tions among the microscopic constituents, energy is not conserved and external 
sources are necessary in order to maintain a non-trivial stationary state: time re- 
versal invarianee is broken and consequently, properties such as the Equilibrium 
Fluctuation-Dissipation relation (EFDR) do not hold. In recent years, a rather 
complete theory for fluctuations in granular systems has been developed for the 
dilute limit, in good agreement with numerical simulations |631I64|. However, a 
general understanding of dense cases is still lacking. A common approach is the 
so-called Enskog correction [63l |65] , which renormalizes the collision frequency 
to take into account the breakdown of Molecular Chaos due to high density. In 
cooling regimes, the Enskog theory may describe strong non-equilibrium effects, 
due to the explicit cooling time-dependence |66| . However it cannot describe 
dynamical effects in stationary regimes, such as violations of the Einstein rela- 
tion [67l |68] • Indeed, as discussed before, violations of Einstein relation comes 
from having neglected coupled degrees of freedom which are expected to be inde- 
pendent at equilibrium. The Enskog approximation is not sufficient to describe 
the effect of such a coupling, because it does not break factorization of velocities: 

N 
p{xi,...xn,vi,...vn) = Px{xi,...,XN)Y\_Pli'^i) (80) 

with pi{v) the single-particle velocity distribution. Such an approximation, 
which predicts the validity of Einstein relation, is revealed to be wrong already 
at moderate densities in granular fluids, as discussed in the following. 

6.1 Model 

A good benchmark for new non-equilibrium theories is provided by the dynam- 
ics of a massive tracer interacting with a gas of smaller granular particles, both 
coupled to an external bath. In particular, taking as a reference point the dilute 
limit, where the system has a closed analytical description |69| , it is shown that 
more dense configurations are well described by a Generalized Langevin Equa- 
tion (GLE) with an exponential memory kernel of the form given in Eq. (j62[) . 
at least as a first approximation capable of describing violations of EFDR and 
other non-equilibrium properties |55j . Here, the main features are: 

• the decay of correlation and response functions is not a simple exponential 
and shows backscattering [70l [Tl] ; 

• the EFDR [721 [73] of the first and second kind do not hold. 

In the model described here detailed balance is not satisfied in general, non- 
equilibrium effects can be taken into account and the correct behavior of cor- 
relation and response functions is predicted. Furthermore, the model has a 
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remarkable property: it can be mapped onto the two-variable Markovian sys- 
tem discussed in the previous section, i.e. two coupled Langevin equations, as 
in Eqs. (p5)) . The dilute limit is then naturally recovered by putting to zero 
the coupling constant between the original variable and the auxiliary one. The 
auxiliary variable can be identified in the local velocity field spontaneously ap- 
pearing in the surrounding fluid. This allows us to measure the fluctuating 
entropy production (see Ea. (|57)) [74], and fairly verify the Fluctuation Rela- 
tion (|5T|) [751 inil [75], a remarkable result, if considered the interest of the 
community [73] and compared with unsuccessful past attempts [77] [7S] . 

The model considered here [SS] is the following: an "intruder" disc of mass 
rriQ = M and radius R, moving in a gas of N granular discs with mass rrii = m 
{i > 0) and radius r, in a two dimensional box of area A = L^. We denote by 
n = N/A the number density of the gas and by 4> the occupied volume fraction, 
i.e. (p = Tr{Nr'^ + R'^)/A and we denote by V (or Vq) and v (or Vi with i > 0) the 
velocity vector of the tracer and of the gas particles, respectively. Interactions 
among the particles arc hard-core binary instantaneous inelastic collisions, such 
that particle i, after a collision with particle j, conies out with a velocity 

v', = Vi - (1 + a) ^ — [{v, - Vj) ■ n]h (81) 

rrii + rrij 

where n is the unit vector joining the particles' centers of mass and a G [0, 1] 
is the restitution coefficient (a = 1 is the clastic case). The mean free path 
of the intruder is proportional to Iq = !/((?' + R)n) and we denote by Tc its 
mean collision time. Two kinetic temperatures can be introduced for the two 
species: the gas granular temperature Tg = m{v'^)/2 and the tracer temperature 
Ttr = M{V^)/2. 

In order to maintain a granular medium in a fluidized state, an external 
energy source is coupled to each particle in the form of a thermal bath [7^[5ni[5T| 
(hereafter, exploiting isotropy, we consider only one component of the velocities): 

m^v,{t) = ^JbV^{t) + f,{t) + Ut)- (82) 

Here fi (i) is the force taking into account the collisions of particle i with other 
particles, and £,b{t) is a white noise (different for all particles), with {^b{t)) = 
and {^b{t)ib{t')) = 2TbjbS{t — t'). The effect of the external energy source 
balances the energy lost in the collisions and a stationary state is attained with 
m,{v1) < Tb . 

At low packing fractions, < 0.1, and in the large mass limit, m/M ^ 1, 
using the Enskog approximation it has been shown [55] that the dynamics of 
the intruder is described by a linear Langevin equation: 

MV = -TeV + £e, (83) 

where £e is a white noise such that 



(£,(i)£,(t'))-2 



h \ 2 ^ 



IbTb + 7, ^^T, 



5,j5{t-t') (84) 
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and 



76 + 7f, with j^ = !iI^-L^^2^^^{l + a) (85) 



where g2{'f+R) is the pan- correlation function for a gas particle and the intruder 
at contact. In this limit the velocity autocorrelation function shows a simple ex- 
ponential decay, with characteristic time M/Te- Time-reversal and the EFDR, 
which are very weakly modified for uniform dilute granular gases |82| 1671 I83j , 
become perfectly satisfied for a massive intruder. The temperature of the tracer 
is computed as T^ = {-ybTb + lg^^^Tg)/rE- For a general study of a Langcvin 
equation with "two temperatures" but a single time scale (which is always at 
equilibrium), see also j84) . 

6.2 Dense case: double Langevin with two temperatures 

As the packing fraction is increased, the Enskog approximation is less and less 
effective in predicting memory effects and the dynamical properties. In par- 
ticular, velocity autocorrelation C{t) = {V {t)V {0)) / {V"^) and linear response 
function R{t) ~ SV{t)/6V{0) (i.e. the mean response at time t to an impulsive 
perturbation applied at time 0) show an exponential decay modulated in ampli- 
tude by oscillating functions [71]. Moreover violations of the EFDR C{t) = R{t) 
(Einstein relation) are observed for a < 1 [68l |41] . 

Molecular dynamics simulations [55j of the system have been performed, 
giving access to C{t) and R{t), for several different values of the parameters a 
and (/), see Fig. S] 

Notice that the Enskog approximation |63l [69] cannot predict the observed 
functional forms, because it only modifies by a constant factor the collision 
frequency. In order to describe the full phenomenology, a model with more than 
one characteristic time is needed. The proposed model is a Langevin equation 
with a single exponential memory kernel [5SJ [SS] [SS] , as in Eq. ([5^ 



MV{t) = - / dt' r(t - t')V{t') + £'{t), (86) 

^— oo 

in this case, 

r(i)= 270(5(0 +7i/rie-*/"i (87) 

and £'{t) = 8o{t)+£i{t), with 



and {£i{t)£o{t')) = 0. In the limit a — > 1, the parameter Ti is meant to tend 
to To in order to fulfill the EFDR of the 2nd kind {£'{t)£'{t')) = ToT{t - t'). 
Within this model the dilute case is recovered if 71 — > 0. In this limit, the 
parameters 70 and To coincide with F^; and T^ of the Enskog theory |69j . The 
exponential form of the memory kernel can be justified within the mode-coupling 
approximation scheme, see [55] for details. 
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Figure 4: (Color online). Semi-log plot of C{t) (symbols) for different values 
oi (f> = 0.01, 0.1, 0.2, 0.33 at a = 0.6. Times are resealed by the mean eoUision 
time Tc. The other parameters are fixed: N = 2500, m = 1, M = 25, r = 0.005, 
R — 0.025, Tb — 1, jb = 200. Times are resealed by the mean collision times 
Tc, as measured in the different cases. Numerical data are averaged over ^ 10^ 
realizations. Continuous lines are the best fits obtained with Eq. ([59]) . Inset: 
C{t) and the best fit in linear scale for </> = 0.33 and a = 0.6. 



The model ^ predicts C = fc{t) and R = fuit) with 



fc{R) 



[cos(wi) + ac{n) sin(a;i)]. 



(89) 



The variables g, ui, ac and an are known algebraic functions of 70, To, 71, ti and 
Tj. In particular, the ratio ac/an = [To - 0(Ti - To)]/[To + 0(ri - To)], with 
fl = ^1 / [{"fo + ji){'yQ / Mti — 1)]. Hence, in the elastic (Ti — )• Tq) as well as in the 
dilute limit (71 —>■ 0), one gets ac = an and recovers the EFDR C{t) = R{t). In 
Fig.|4]the continuous lines show the result of the best fits obtained using Eq. (|89p 
for the correlation function, at restitution coefficient a = 0.6 and for different 
values of the packing fraction 0. The functional form fits very well the numerical 
data. A fit of measured C and R against Eqs. ([M]) . together with a measure of 
(V^), yields five independent equations to determine the five parameters entering 
the model. We used the external parameters mentioned before, changing a or 
the box area A (to change (jj) or the intruder's radius R, in order to change 
keeping or not keeping constant 7^ ^ 1/Iq — >■ (indeed different dilute limits can 
be obtained, where collisions matter or not). Fits from numerical simulations 
suggest the following identification for the parameters: 70 '^ Te, Tq ^ Ttr and 
Ti ~ Tfc. The coupling time ri increases with the packing fraction and, weakly, 
with the inelasticity. In the most dense cases it appears that 71 ~ 7^ ex <p: such 
an observation however docs not hold in the dilute limit at constant collision 
rate, where 71 ^^ ^ 7^". It is also interesting to notice that at high density 
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Ttr ^ Tg ^ T^ , which is probably due to the stronger correlations among 
particles. Finally we notice that, at large 0, Ttr > T^, which is coherent with 
the idea that correlated collisions dissipate less energy. 

Looking for an insight of the relevant physical mechanisms underlying such 
a phenomenology and in order to make clear the meaning of the parameters, it 
is useful to map Eq. (j86p onto a Markovian equivalent model by introducing an 
auxiliary field (see |53l 155] for details on how the mapping is achieved): 



V 



MV = -TE{V-U) + ^/2r]^g& 
M'U = -T'U - TeV + ^/2T'Tb8u, (90) 

where £y and En are indipendent white noises of unitary variance and we have 
exploited the numerical observations discussed above (i.e. Tq ^ Tg , Ti ^ Tf, and 

2 2 

70 ~ Tb), while F' = — and M' — ^iill, Jn the chosen form (|90|) . the dynamics 
of the tracer is remarkably simple: indeed V follows a memory less Langevin 
equation in a Lagrangian frame with respect to a local field U, which is the local 
average velocity field of the gas particles colliding with the tracer. Extrapolating 
such an identification to higher densities, we are able to understand the value 
for most of the parameters of the model: the self drag coefficient of the intruder 
in principle is not affected by the change of reference to the Lagrangian frame, 
so that 7o ~ Te', for the same reason Tq ^ Ttr is roughly the temperature of 
the tracer; ti is the main relaxation time of the average velocity field U around 
the Brownian particle; 71 is the intensity of coupling felt by the surrounding 
particles after collisions with the intruder; finally Ti is the "temperature" of the 
local field U, easily identified with the bath temperature Ti ^ Ti,: indeed, thanks 
to momentum conservation, inelasticity does not affect the average velocity of 
a group of particles which almost only collide with themselves. 

6.3 Generalized FDR and entropy production 



The model discussed above, Eq. (|90)) . is able to reproduce the violations of 
EFDR, as show in Fig. [Sj which depicts correlation and response functions 
in a dense case (elastic and inelastic). In the inelastic case, deviations from 
EFDR R(t) = C(t) are observed. In the inset of Fig. [5] the ratio R{t)/C{t) 
is also reported. As shown in Section ^ a relation between the response and 
correlations measured in the unperturbed system still exists, but - in the non- 
equilibrium case - must take into account the contribution of the cross correla- 
tion {V{t)U{0)), i.e.: 

Rit) = aC{t) + b{V{t)U{0)) (91) 

with a = [1 — ji/M{Tq — Ti)Qa] and b = (Tq — Ti)Qb, where fla and Q^ are 
known functions of the parameters. At equilibrium, where Tq = Ti, the EFDR 
is recovered. 

An important independent assessment of the effectiveness of model in Eq. (|90|) 
comes from the study of the fluctuating entropy production |74j which quantifies 
the deviation from detailed balance in a trajectory. Given the trajectory in the 
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Figure 5: (Color online). Left: correlation function C{t) (black circles) and 
response function R{t) (red squares) for a = 1 and a = 0.6, at = 0.33. 
Continuous lines show the best fits curves obtained with Eqs. (|89)) . Inset: the 
ratio R{t)/C{t) is reported in the same cases. Right: Check of the fluctuation 
relation (|93| in the system with a = 0.6 and </> = 0.33. Inset: collapse of the 
rescaled probability distributions of St at large times onto the large deviation 
function. 



time interval [0,i], {^(s)}q, and its time-reversed {IV{s)}q = {—V{t- 



s)y 



0' 



as 



discussed in section 
takes the form \M\ 



shown that the entropy production for the model (|M|) 



log 






7o 



1 



1 



ds V{s)U{s) 



(92) 



This functional vanishes exactly in the elastic case, a = 1, where equipartition 
holds, Ti = To, and is zero on average in the dilute limit, where (VU) = 0. 

Formula (|92p reveals that the leading source of entropy production is the 
energy transferred by the "force" jqU on the tracer, weighed by the difference 
between the inverse temperatures of the two "thermostats" . Therefore, to mea- 
sure entropy production, we need to measure the fluctuations of U: from the 
above discussion we have for U the interpretation of a spontaneous local velocity 
field interacting with the tracer. Therefore it can be measured doing a local av- 
erage of particles' velocity in a circle of radius / centered on the tracer. Details 
on how to choose in a reliable way the proper I are given in |55| . for instance 
following such procedure, in the case </> = 0.33 and a = 0.6, we estimate for the 
correlation length / '^ 9r ~ GIq- Then, measuring the entropy production from 
Eq. (inn along many trajectories of length t, we can compute the probability 
P(Et = x) and compare it to P{T,t ~ —x), in order to verify the Fluctuation 
Relation 

P{^t = x) 



log- 



(93) 



Pi^t = -x) 

In the right frame of Fig. [5] we report our numerical results. The main frame 
confirms that at large times the Fluctuation Relation (|93p is well verified within 
the statistical errors. The inset shows the collapse, for large times, of log P{Yit)/t 
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onto the Cramer function S{y = ^), introduced in Section |31 Notice also 
that formula ((92|) does not contain further parameters but the ones already 
determined by correlation and response measure, i.e. the slope of the graph is 
not adjusted by further fits. Indeed a wrong evaluation of the weighing factor 
(1/To - l/Ti) « (l/Ttr - 1/Tb) or of the "energy injection rate" joU(t)V{t) in 
Eq. ([^ could produce a completely different slope in Fig. [S] (right frame). 

To conclude this section, we stress that velocity correlations {V{t)U{t')) be- 
tween the intruder and the surrounding velocity field are responsible for both 
the violations of the EFDR and the appearance of a non-zero entropy produc- 
tion, provided that the two fields arc at different temperatures. We also mention 
that larger violations of EFDR can be observed using an intruder with a mass 
equal or similar to that of other particles [68], with the important difference 
that in such a case a simple "Langevin-like" model for the intruder's dynamics 
is not available. 



7 Conclusions and perspectives 

We have reviewed a series of recent results on fluctuation-dissipation relations 
for out-of-equilibrium systems. The leitmotiv of our discussion is the impor- 
tance of correlations for non-equilibrium response, much more relevant than in 
equilibrium systems: indeed the Generalized Fluctuation-Dissipation relation 
discussed in Section |31 El [HI deviates from its equilibrium counterpart for the ap- 
pearance of additional contributions coming from correlated degrees of freedom. 
This is the case, for instance, of the linear response in the diffusion model on a 
comb-lattice, where - in the presence of a net drift - the linear response takes a 
non- negligible additive contribution, see Eq. ([FT]) . The same occurs in the gen- 
eral two-variable Langevin model: additional contribution to the equilibrium 
linear response appears when the main field V is coupled to the second field C/, 
which only appears out-of-equilibrium, see Eq. (j7ip . The granular intruder is 
a realistic many-body instance of such a coupling scenario, where we have seen 
the strong influence of correlated degrees of freedom to the linear response, see 
Eq. 0. 

Remarkably, in some cases, one may explicitly verify that the coupled field 
which contributes to the linear response in non-equilibrium setups is also in- 
volved in the violation of detailed balance: such a violation is measured by 
the fluctuating entropy production. The connection between non-equilibrium 
couplings and entropy production has been discussed for general coupled linear 
Langevin models, see Eq. ([7^ and then verified for the entropy production of 
the granular intruder, Eq. ([92)) . 

In conclusion many results point in the same direction, suggesting a gen- 
eral framework for linear-response in systems with non-zero entropy production. 
Even in out-of-equilibrium configurations, a clear connection between response 
and correlation in the unperturbed systems exists. A further step is looking 
for more accessible observablcs which could make easier the prediction of linear 
response: indeed, both proposed formula, Eq. ([T3|) and Eq. (pO| . require the 
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measurement of variables which, in general, depend on full phase-space (micro- 
scopic) observables and can be strictly mo del- dependent. Such a difficulty also 
explains why, in a particular class of slowly relaxing systems with several well 
separated time-scales, such as spin or structural glasses in the aging dynam- 
ics, i.e. after a sudden quench below some dynamical transition temperature, 
approaches involving "effective temperatures" have been used in a more satis- 
factory way [87j . More recent interpretations of the additional non-equilibrium 
contributions has been proposed in [55] , but its predictive power is not yet fully 
investigated and represents an interesting line of ongoing research. 
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